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PALM DISTRIBUTIONS OF WAVE CHARACTERISTICS 
IN ENCOUNTERING SEAS 
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Lund University, Lund University and University of North Carolina 

Distributions of wave characteristics of ocean waves, such as wave 
slope, waveheight or wavelength, are an important tool in a variety 
of oceanographic applications such as safety of ocean structures or 
in the study of ship stability, as will be the focus in this paper. We 
derive Palm distributions of several wave characteristics that can be 
related to steepness of waves for two different cases, namely for waves 
observed along a line at a fixed time point and for waves encountering 
a ship sailing on the ocean. The relation between the distributions 
obtained in the two cases is also given physical interpretation in terms 
of a "Doppler shift" that is related to the velocity of the ship and the 
velocities of the individual waves. 

1. Introduction. The study of wave characteristics such as crest height, 
wavelength and wave slope is important in various ocean engineering ap- 
plications. Examples can be found in the design of ocean structures, such 
as oil platforms, or in the design of sea walls that should prevent cities, or 
sometimes whole countries, from flooding. Another important application, 
motivating this paper, is that of ship stability. Because high and steep waves 
encountering a ship may cause structural damage and even capsize smaller 
vessels, it is important to know the distribution of, for example, wavelength, 
waveheight and wave slope of such waves. Based on these distributions, cap- 
size probabilities can be computed and used as a risk measure for existing 
vessels or as a tool in the design of new ones. 

A sea surface can be seen as a sequence of apparent waves. By an apparent 
wave is meant the part of the sea record between two consecutive upcrossings 
of the still water level and the downcrossing in between these upcrossings is 
called the center of the wave. The purpose of this paper is to derive exact 
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distributions of steepness related wave characteristics of the apparent waves 
for two different cases. The first case, the spatial case, concerns properties 
of the waves observed, in space, along a fixed line at a fixed time point 
and the second case, the encountered case, involves properties of the waves 
encountering a ship on the ocean. More precisely, the characteristics of those 
waves that overtake a ship sailing with constant velocity along a straight 
line will be investigated because such waves are considered as particularly 
dangerous for ship stability. 

Distributions of wave characteristics are defined as Palm distributions. 
Such distributions are intimately related to level-crossings of the underlying 
process; in our case representing the sea surface, and Rice's formula proves 
useful in their evaluation. Although the main result of this paper is the 
derivation of several distributions for the two different cases stated above, 
a nice by-product is that the relation between the two cases can be given 
interpretation in terms of wave velocities. In physical terminology this can be 
expressed as a Doppler shift and is caused by the fact that the sea surface is 
observed by a moving observer — the ship. Due to this interpretation, similar 
results can be expected in many other applications where the distribution 
of moving objects, observed from a moving observer are studied. As an 
example, the distribution of the size of the storms you will meet when sailing 
will be related, by a Doppler shift, to the distribution of storm sizes that 
can be observed from a satellite image. 

The organization of this paper is as follows. First, the Gaussian sea model 
used to evaluate the distributions for the wave characteristics is introduced. 
Then it is shown how distributions defined by Palm distributions can be 
computed by using generalized versions of Rice's formula. After a short 
section on wave velocities, the theory is exemplified by calculation of the 
distribution of wave slope for spatial and encountered waves. Thereafter, a 
more intricate example concerning waveheight and wavelength of apparent 
waves is investigated, and finally, the derived distributions are evaluated 
numerically for a unidirectional Gaussian sea. 

2. Gaussian sea model. Let W(x,t) be the sea surface elevation at lo- 
cation x and time t. Further denote its partial derivative on t by Wt{x,t) 
and its first and second partial derivative on x by W x {x,t) and W xx (x,t), 
respectively. Sometimes, when no misunderstanding can be made, the nota- 
tion W(x) = W(x,0) is used. Throughout this paper, W(x,t) is modeled as 
a zero mean, stationary Gaussian field with a directional spectrum having 
spectral density S(u>, 6), where u is an angular frequency and 9 an angle rep- 
resenting wave direction. These quantities satisfy u > and 8 € [0,2-7r]. For 
sea waves such a parametrization of the spectral density is possible because 
the angular frequency u and the wave number k are related by physical 
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dispersion laws; see [4] for further reference. In particular, for deep water, 

2.1 K=— , 

9 

where g is the gravitational constant. Writing uj = 2n/T and k = 2-k/L, 
where T is the wave period and L the wavelength, the dispersion relation 
can also be expressed in terms of periods and wavelengths, namely 

L = ^-T 2 . 

2vr 

Due to the dispersion relation, we thus have the following relation between 
the covariance function R(£, r) = E[W(a;, t)W(x + £, t + r)] and the spectrum 

pod r2iT / f.2 \ 

(2.2) R(£,r) = J J caB\WT + —ZcoB0jS(u},6)dBdw. 

Due to this relation variances and covariances of the process and its deriva- 
tives can be expressed in terms of the spectral moments defined as 

poo r2iT 

(2.3) Xij= / / (k(uj) cos eyu j S(uj, 6) dOdu. 

Jo Jo 

Note that Var(W(0,0)) = R(0,0) = A 00 , Var(W a: (0,0)) = -%(0,0) = A 20 , 
Vax(W t (0,0)) = -Rtt(0,0) = A 02 and 0^(^(0,0), W t (0,0)) = -R^ T (0,0) = 
An. 

3. Palm distributions and Rice's formula. Let iV|([0,T]) be the number 
of times that a process Z(t) of one parameter takes the value u in the 
interval [0,T], and let N^([0,T],E) be the number of times the process 
takes the value u and at the same time a statement E about the process, 
its derivatives or another process, is satisfied. For example E could be the 
statement u Zt(t) < z," where Zt(t) is the derivative of Z(t). A Palm measure 
is defined by the following ratio of intensities 

r , n p n fFl _ E[iVl([0, !],£)] 

(3 - 1} z{ ) ~ m¥^)\ ' 

provided that the expected number of u-crossings, E[iV^([0, 1])], is finite. If 
the process Z(t) is ergodic the intensities in (3.1) can be computed as sample 
averages and, with probability one, 

P z {E)-\^ NmT]) ■ 

Thus the Palm distribution has the empirical interpretation that it is the 
long-term proportion of it-crossings by the process Z for which E is satisfied. 
Note that a sufficient condition for ergodicity for a Gaussian process is that 
it possesses a spectral density function; see [7], page 157. 
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The intensities of crossings in (3.1) can be computed by using Rice's 
formula, first studied by [8] and [15]. In its simplest form, when Z is a zero- 
mean stationary Gaussian process having almost surely differentiable sample 
paths (realizations), it reads 



tively. There is a vast literature on generalized versions of Rice's formula; see, 
for example, [1, 3, 5, 9, 18]. It turns out that the formulation (3.3) extends 
in a natural way so that the intensity of crossings satisfying a statement 
E can be computed. For example, if Z(t) has almost surely differentiable 
sample paths and the event E satisfies certain regularity conditions, then 
by [9], Lemma 7.5.2. 



of it can be used to compute the Palm distribution (3.1). However, some 
of the Palm distributions that we will consider include crossings by random 
vector fields and therefore slightly more general forms of Rice's formula are 
required. The exact version that will be used is due to [13], Theorem 9.6, 
which is a generalization of a theorem in [3] . 

Remark 3.1. Notation like E[|Z f (t)| | Z(t) = u}f Z (t){u) is used exten- 
sively throughout this paper. This should be interpreted in the sense of (3.3), 
that is, as an integral. In particular, this requires that the joint density of 
Zt(t) and Z(t) exists, so that (Zt(t),Z(t)) has a nondegenerate Gaussian 
distribution. For higher dimensions the notation should be understood in a 
analogous fashion. 

Remark 3.2. It should be pointed out that the Gaussian assumption 
is not necessary for Rice's formula (3.3) to hold, although the conditions for 
its validity in some cases can be quite intricate. However, if one is content 
to use (3.3) for almost every u, rather than for each specific u, these condi- 
tions simplify substantially. In [10] sufficient conditions for applying almost 
everywhere results to fixed levels are given. 



(3.2) 




It can also be written on integral form, namely 
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4. Wave velocities. We will see that encountered distributions of wave 
characteristics are closely related to velocities of the individual apparent 
waves. Before defining velocities for these random waves we start with ve- 
locities of deterministic waves, and in particular, study the implications of 
the dispersion relation for the velocities of such waves. 

Consider a deterministic cosine-wave c(x,t) with angular wave frequency 
u! and wave number k, namely 



where A is an amplitude and 4> a phase. For fixed t this is a wave in space 
and for fixed x, a wave in time. If both t and x are varying it is a wave 
traveling in the direction of the negative x-axis. The (phase) velocity of the 
wave can be expressed as 



Note that the convention here is that a wave travels in the direction of 
the positive x-axis if it has positive velocity and vice versa. Now for this 
deterministic wave, uj and k are linked by the dispersion relation (2.1) and 
consequently 



Thus for ocean waves the dispersion relation implies that longer waves are 
faster than shorter ones and, because the spatial slope of the wave at a zero- 
crossing is proportional to k, it also implies that waves with a high spatial 
slope are slow. 

Velocities for waves in irregular seas can be defined in a similar fashion, 
and has been studied by [2, 4, 12, 14]. In this case, however, one is interested 
in the velocities of the apparent waves rather than velocities of single har- 
monics. Recall that an apparent wave is the part of the sea record between 
two consecutive upcrossings of the still water level, and that the downcross- 
ing in between is termed the center of the wave. Let Xi > be the positions 
of centers of waves observed in W(x,0), that is, along the x-axis at time 
zero. Due to the time variability of the sea surface, the centers change their 
positions with time and at t = the velocity Vi that the ith wave moves 
at can be evaluated by means of Vi = — w^ixfo) > see ^ or exam pl e i [4]- The 
variability of the velocities Vi will be described by the Palm distribution 
defined as follows 



c(x, t) = Acos(ujt + kx + (j)) 




c t (x,t) 
c x (x, t) 




V = 



K 




(4.2) 



F v {v) 



E [number of Xj < 1 such that Vi < v] 



E [number of Xj < 1] 
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For the Gaussian sea this becomes 




V>-^ + a7A 20 



v — v 



) 




where a 2 = A02 — and v is the average wave velocity; see [14]. 

Remark 4.1. The velocity defined by V(x,y,t) = - wlfx^t) » ( x >v) G 
M 2 , can be seen as the local velocity field of a random, moving, surface. This 
velocity is defined at any point (x,y) € M 2 at any time t. For a thorough 
study of different velocity concepts; see [4]. 

5. Explicit evaluation of Palm distributions. Earlier it was shown how 
Palm distributions can be expressed by Rice's formula, given by (3.3) and 
(3.4). These formulas involve computation of multivariate normal expecta- 
tions. Sometimes these have to be evaluated numerically, but in some cases 
explicit forms can be given. One of these explicit cases will be investigated 
next in some detail. 

In the following lemma, a stationary Gaussian vector valued process (Z(t), 
Zt(t), Y(t)) will be considered. In order to simplify notation, introduce the 
random variables Y u and Zf having joint distribution function iV«,z u (^^) 
defined by the measure (3.1), namely 



(5.1) F Y u >Z u( y , z) = P{Y U < y, Z? <z) = Pg(Y(t) < y, Z t (t) < z). 



These random variables are a special case of the Slepian model process; see 
[11] for a detailed presentation, and one may think of them as Y(t) and 
Zt(t) observed at a randomly chosen -u-crossing of the process Z(t). The 
next lemma gives an explicit representation of these variables. 

Lemma 5.1. Suppose that (Z(t),Zt(t),Y(t)) is a stationary Gaussian 
vector valued process such that the density of (Z(t), Z t (t),Y(t)) exists. As- 
sume further that the sample paths of Y(t) and Z(t) are a.s. continuous 
and a.s. continuously differentiable, respectively. Let mz = E(Z(t)), my = 
E(y(<)) and let the covariance matrix E of the vector (Z(t), Zt(t),Y(t)) be 
given by 



Then the variables Y u ,Zf, with distribution defined by (5.1), have the fol- 
lowing distributional representations: 



(5.2) 
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where U is standard normal and R is independent of U , with a double 
Rayleigh distribution, that is, probability density f(r) = ^e _r2 / 2 . Here m u = 
m Y + ^Pzy(u — m z ) and p ZY , Pz t Y are correlations between Y(0) and 
Z(0), Zt(0), respectively. 

Proof. The proof is easy using the definition of Palm probability, Rice's 
formula (3.2) and (3.4) and the Gaussian regression formulas. □ 

This lemma thus shows that the derivative Zt at crossings has a Rayleigh 
distribution as opposed to the derivative at any point, where it has a Gaus- 
sian distribution. Furthermore, the distribution of another Gaussian process 
Y, correlated with Z(t), can be represented as the sum of a Rayleigh and 
Gaussian variable at the crossing point. The corresponding distributions, but 
only observed at up- or downcrossings of Z(t), are easily obtained by using 
the representation from Lemma 5.1. These distributions are summarized in 
the following lemma. 

Lemma 5.2. Let Zf \ Zf > denote the random variable obtained by 
conditioning Z? on Zf > and let Y u \Z?>0, Z? \ Zf < and Y u \Z?<0 
be defined analogously. Using the same notation as in Lemma 5.1, we then 
have the following distributional representations: 

Z?\Z?>0 = a Zt R+, 

Y u \Z?>0 = m u + a Y (p ZtY R + + yl - p\ Y - p\ Y U) 

and 

Z?\Z?<0 = a Zt R-, 

Y u \Z?<0 = m u + o-y( PZtY R- + yjl - p 2 ZY - p% tY U), 

where R + and R~ , both independent of U , have densities fn+(r) = re _r2//2 , 
r > 0, and (r) = — re~ r I 2 , r < 0, respectively. 

This result will be used to obtain explicit representations of the distribu- 
tions of the slope observed at wave centers and encountered wave centers. 
When evaluating such distributions, the following lemma that can be proved 
by straightforward calculations is useful. 

Lemma 5.3. Let U and R + be independent standard Gaussian and Rayleigh 
distributed variables. Then, for a > and any x, 

P(aU + bR+> X ) = $ (-*) + e^l^ ■ U ( X - h X 

where a = \J a 2 + b 2 and $ is the distribution function of a standard Gaus- 
sian variable. 
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6. Distribution of slope. Consider the sea surface W(x, t) . In this section 
the distribution of the slope W x (x,t) will be computed, but only for values 
of x and t chosen in a careful manner. The cases considered are: 

1. Spatial case. Choose x to be points where the process W(x,0) has a 
downcrossing of the zero level. This means that only values of x satisfying 
W(x,0) = and W x (x,0) < are chosen. In other words, x are locations 
of the wave centers in W(-,0). These values of x and t = will lead to a 
Palm distribution of wave slopes observed at centers of waves. 

2. Encountered case. Consider a ship sailing with constant velocity v on the 
ocean. If it at time zero is at position x = 0, then the sea elevation at 
the center of gravity at time t is W(vt,t). Choose t so that W(vt,t) = 0, 
W x (vt,t) < and -j^W(vt,t) < 0, then t can be interpreted as a time 
when the ship is overtaken by a wave center. Values of t chosen in this 
manner, and x = vt, lead to a Palm distribution of wave slopes observed 
at overtaking centers of waves. 

Remark 6.1. By slope is here meant the water surface wave slope. This 
should not be confused with the term velocity slope sometimes used in the 
marine sciences. 

6.1. Spatial distribution of slope. To derive the distribution of the slope 
observed at centers of waves, we start by defining it by means of a Palm 
distribution. Therefore, let Xj > be the locations of the centers of waves, 
that is, locations of the zero-downcrossings, in W(x,0). Each of these wave 
centers can be associated with a slope W x (xi,0) and the following Palm 
distribution for the slope at centers of waves, denoted by F^ ace (ui), can be 
defined as 



F space ( 



W ) 

^ E[number of X{ < 1 such that W x (xi,0) < w] 

E [number of X{ < 1] 

= E[#{x£ [0,1]; W(x,0)=0, W x (x,0) < 0,W x (x, 0)<w}] 
E[#{x e [0, l];W(x, 0) = 0, W x (x, 0) < 0}] 

where the second equality holds because point of downcrossing of 

W(x,0). Note that by ergodicity this distribution can be interpreted as 

,-,sDace/ \ , number of Xi < X such that W x (xi,0) < w 

F w ( w ) = lim Z 1 v > 

Wx x^oo number of Xi < X 

that is, as the proportion of wave centers with associated slope less or equal 
to w in an infinitely long realization of the process. 

The following theorem gives an expression for the Palm distribution in 
this case, and also states how it can be evaluated. 
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Theorem 6.1. IfW(x,0) is a stationary Gaussian process having a.s. 
continuously differentiate sample paths then 



Proof. The first statement of the proof follows by applying the gener- 
alized Rice's formula (3.4) to the numerator and denominator of (6.1). For 
the numerator, (3.4) is used with E = "W x (x, 0) < 0, W x (x, 0) < w" and for 
the denominator E = u W x (x,0) < 0." 

To prove (6.3), Lemma 5.1 is used with Z = W(x, 0). Then Z x = W x (x, 0), 
and the Slepian variable Z® is the derivative W x (x,0) observed at zero- 
crossings of the process Z = W(x, 0). Using the variable Z®, the distribution 
for the slope W x (x,0) observed at zero-downcrossings can be expressed as 



that is, as the distribution of the slope observed at zero-crossings conditional 
on a negative slope at the crossing. Now (6.3) follows by Lemma 5.2. □ 

6.2. Encountered distribution of slope. Consider a vessel sailing along 
the x-axis in the positive direction with constant speed v, having its center of 
gravity at time zero at x = 0. Disregarding the interaction between the ship 
and the waves, the sea elevation at the center of gravity of the ship is given 
by Z(t) = W(vt,t). In the following, we will refer to Z(t) as the encountered 
sea. Now let ti > be the times when the vessel is overtaken by a wave 
center. Characteristic for such times is that the center of gravity is passing 
through the still water level, Z{ti) = 0, the slope of the wave is negative, 
W x (vti,ti) < and the encountered sea Z(t) has an upcrossing Z t {ti) > 0, 
meaning that the center of wave is overtaking the ship. Because each ti 
can be associated with the slope W x (vti,ti) of the iih overtaking wave, the 
following Palm distribution of the slope at overtaking waves, F^(w), say, 
can be formed 




E[W g (0,0)-l {Wc( o i0) < 1l , } | ^(0,0) = 0]/hW°) 
E[^(0,0)-|^(0,0) = 0]/ W(0i0) (0) 



where x = max(— x,0). Moreover 
(6.3) F^ cc (w) = 




where R is a random variable having density f R - (r) 



re r2 / 2 , r < 0. 



F^ cc (w)=P(Z° x <w\Z° x <0) 




(6.4) 



E[number of ti < 1 such that W x (vti,ti) < w] 



E [number of ti < 1] 

E[#{t g [0, 1]; Z(t) = 0, Z t (t) > 0, W x {vt, t) < 0, W x (vt, t) < w}} 
E[#{t G [0, 1]; Z(t) = 0, Z t (t) > 0, W x {vt, t) < 0}] 
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where the second equality is due to the definition of ij. As before, the Palm 
distribution can be interpreted in the following frequency fashion: 

,_,„„„. . .. number of ti < T such that W x (vti,ti) < w 

FyTr (w) = hm — . 

w * v 1 t^oo number of U < T 

Remark 6.2. Despite the similarities of expressions (6.1) and (6.4), 
there is an important difference arising from the fact that crossings from 
two different processes are counted, namely crossings from W(x, 0) in the 
spatial case and from W(vt,t) in the encountered case. From a statistical 
point of view this means sampling from two different populations. 

The following theorem, concerning the Palm distribution of slope for en- 
countered waves, is an analogue to Theorem 6.1. 

Theorem 6.2. Let W(x,t) be a stationary Gaussian process having a.s. 
continuously differ -entiable sample paths. Then Fyy c (w) can be expressed as 

(6 51 F enc ( ) = E [ Zf (°) +1 {W J; (o,o)<o,^(o,o)<^} I Z{0) = 0]/z(o)(0) 
{ '> WAW> E[Z t (0)+l {Wxm<0} \Z(0) = 0]f m (0) ' 

where x + = max(x, 0) . Moreover, Fyy c (w) = 1 if w > and for w < 
= ( $ i n W 5 ) 

l~pZ t Y\ \J\oJl-ol , A J 



(6.6) 



'A 2 o(l-p| t y 



PZ t YW 



'^AsoU-p^y)' 

where $ is the distribution function of a standard normal variable and 

v\ 2 o + An 



(6-7) p Zt Y 



v /A 2 o('w 2 A2o + 2v\n + A02) 



Proof. Expression (6.5) follows, in a similar fashion as in the proof 
of Theorem 6.1, by applying the generalized Rice's formula (3.4) to the 
numerator and denominator in (6.4). 

To prove the second statement, identify the process Z(t) of Lemma 5.1 
with the encountered process W(vt,t) and let Y(t) = W x (vt,t). With this 
choice of variables the covariance matrix, expressed in terms of spectral 
moments defined by (2.3), becomes 



(6.8) S = ( v 2 X 20 + 2v\ n + A 02 

vX 2 o + An 
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In this setting Z 4 ° is the time derivative of the encountered process observed 
at its zero-crossings, and Y° is the spatial slope observed at zero-crossings of 
the encountered process. The distribution in (6.5), however, is the distribu- 
tion of the spatial slope, but only at upcrossings of Z(t) such that the spatial 
slope is negative. Using the variables Z® and Y°0 this can be expressed as 

F^{w) = P{Y° <w\Z?>0,Y°< 0) 

P(Y°<w\Z$ > 0) 

P(Y° < 0|Z? > 0) ' W<0, 
1, w>0. 

Use of Lemma 5.2 now gives the following representation 



P{^(p ZtY R + +yjl- p 2 ZtY U) < W) 

i , = , w < 0, 

*W*W~\ P(^- ( pZtY R+ + ^l-pl Y U)<0) 

1, ^>o, 

with pz t Y given by (6.7). Thus the slope observed at an encountered over- 
taking center of wave can be evaluated by computing probabilities for a sum 
of a Rayleigh and a Gaussian random variable. Thus, (6.6) holds by Lemma 
5.3. □ 

The interpretation of (6.5) is that it is the distribution of the slope of 
the waves that overtake a ship sailing on the ocean. If a wave overtake the 
ship, then it must necessarily be traveling at a higher speed than the ship. 
Thus intuition suggests that in some way the distribution (6.5) should be 
related to the relative velocities of the individual waves and the ship. The 
exact relation is given by the following corollary, enabling us to interpret 
the encountered distribution in terms of physical quantities. 

Corollary 6.1. Let V = —W t (0, 0)/W x (0,0) be the local velocity at 
(x,t) = (0,0) and let v denote the velocity of the ship. Then 

(6,, ^ W 

E [(V - v)+W x (0, 0)-1 { ^(Q,0)<^} I W(0,0) = 0]JV(o,o)(0) 
E[(V-v)+W x (0,0)-\W(0,0) = 0]f wm (0) 

Proof. The assertion of the corollary follows easily from (6.5) due to 
the following equality 

Z t (0) + l {Wx{m<0} = (vW x (0, 0) + W t (0, 0)) + l { ^(o,o)<o} 
= (W x (0,0)(v-V)) + l {Wx{0 , 0)<0} 

(6.10) 

= W x (0,0)-(v-V)- 



W x (0,0)-(V-v 



,+ 
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and the fact that Z(0) = W(0, 0). □ 

The encountered distribution given by (6.5) or (6.9) differs from the cor- 
responding spatial distribution (6.2) only by the term (V — v) + that enters 
into the expectation. Thus the spatial measure for slope can be transformed 
to the encountered measure for slope by entering the velocity term (V — v ) + . 
The physical interpretation of this is that in the encountered approach not 
every center of wave is sampled but only those with velocity greater than 
the ship. The velocity factor is thus a consequence of the fact that the sea 
surface is observed from a moving observer, that is, the ship. In physics, such 
a phenomenon would be termed a Doppler shift. In the following, we will 
see that the Doppler shift transformation of the distribution is not limited 
to the distribution of the spatial slope. In fact, by the way the measures 
are defined, it extends to all wave characteristics of overtaking encountered 
waves. 

7. Distribution of waveheight and half-wavelength. Let x\ be the loca- 
tion of a wave center, that is, a downcrossing of the level zero, in W(x,t) 
for some time t. Then x\ is characterized by the fact that W(x\,t) = and 
W x (xi,t) < 0. Further, let X2 be the distance from the downcrossing x\ to 
the closest local maximum before x\ so that W x (x\ — X2,t) = 0, W xx (x\ — 
X2,t) < and W x {x\ — s,t) < 0, Vs £ (O,^)- The latter condition guaran- 
tees that this really is the closest maximum before x±. Similarly define X3 
to be the distance from x\ to the closest local minimum after x\. Hence 
W x (x! +a? 3 ,t) = 0, Wxx(xi +x 3 ,t) > and W x {x\ - s,t) < 0,Vs G (0,aj 3 ). 
The definitions of X2 and X3 are shown in Figure 1. In the same figure, it is 
also indicated that the notation H2 = W(x\ — X2, t) and H3 = W(x± + £3, t) 
will be used. 

The task is to obtain an expression for the joint Palm distribution of the 
distances from the wave center to the two closest local extrema and their 
heights (x2,xs, H2,H^). As before two cases, corresponding to choosing wave 
centers in two different ways, are considered. The first is a spatial case where 
the Palm distribution is derived for wave centers in W(x,0) and the second 
is an encountered case where the corresponding distribution is derived for 
wave centers that overtake a ship. 

By dividing the waveheight by the half- wavelength a measure of the steep- 
ness at the wave center is obtained, hence motivating the study of these 
quantities from a ship stability perspective. 

7.1. Spatial distribution of heights and distances. In this section, an ex- 
pression for the joint Palm distribution of (x2,x^,H2,H^) at downcrossings 
in space is derived. To do so, we first define the appropriate Palm distribution 
and then express it in terms of level crossings of the process W(x) = W(x,0). 
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Let (x\)i > be positions of downcrossings in W{x). With each down- 
crossing associate distances (x2,x 3 )i, say, to the closest local extrema and 
form the following Palm distribution: 

7-ispace / \ 

F xlx 3 ,H 2 ,H 3 ( r ' S > U > W ) 

E[number of {x\)i < 1 such that (X2, £3, H2, H 3 )i < (r,s,u,w)] 

E [number of < 1] 

In order to express the Palm distribution in terms of level crossings of the 
process W, note that for each crossing (xi)j satisfying (x2, £3, H2, H 3 )i < 
(r, s, u, w) the triple (x±,X2, X3) is characterized by the following properties: 

(Al) (xi,X2,xs) E B, where the set B is defined by 

B = {x G K 3 ; < x x < 1, < x 2 < r, < x 3 < s}, 

(A2) W( Xl ) = W x ( Xl - x 2 ) = W x (x! + x 3 ) = 0, 

W xx (xi - x 2 ) < 0, 
(A3) W XX { X1 + x 3 ) > 0, 

W x (s) < Vs £ T x , where r x = {s 6 K; si - S2 < s < 11 + x 3 }, 

(A4) VF(xi-x 2 )<n, W(xi + x 3 ) < io. 




Fig. 1. Definition of waveheights and distances. If x\ is the location of a downcrossing 
in W(x,i) for fixed t, then X2 is defined to be the distance to the first local maximum before 
Xi and similarly X3 is the distance to the first local minimum after xi. H2 is defined as 
the sea elevation W(xi — X2,t) and H3 as the sea elevation W(xi + X3,t). 
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Defining 

(7.1) C(x)=Z(x 1 ,x 2 ,x 3 ) = (W(x 1 ),W x (x 1 -x 2 ),W x (x 1 + x 3 )) 
it thus holds that 

7-ispace / \ 

, « F x 2 ,x 3 ,H 2 ,H 3 (r,s,u,w) 

7.2 

_ E[#{x £ B;£(x) = 0, (A3) and (A4) satisfied}] 
E[#{xx G [0,l];W( Xl ) = 0,W x ( Xl ) <0}] 

The formulation (7.2) is suitable when it comes to computation of the ex- 
pectations by means of Rice's formula. In particular, the following theorem 
states how the expectation in the numerator can be computed. 



Theorem 7.1. Let W:M— »R be a stationary Gaussian process having 
a.s. twice continuously differentiable sample paths. Assume that the spec- 
trum of W has a continuous component and that the variance of the number 
of zeros W x (y) = 0, y E [0,1] is finite. From W define a process £ by the 
relation (7.1) and a set B as in (Al). Furthermore, define a vector valued 
process Y = (W, W x , W xx ), a set T x as in (A3) and let g(Y, x) be an indicator 
function defined by 

g(Y,x) = l{w x (s)<oys&r x }l{w{x 1 -x 2 )<u}^{w{x 1 +x i )<w} 

X 1 {W xx {x 1 -x 2 )<0} 1 {W xx {x 1 +x 3 )>0}- 

Then, writing Ny(B,g) = #{x £ B;£(x) = 0,g(Y,x) = 1} and det(£ x (x)) for 
the Jacobian determinant of£(x), 

E[N^(B,g)} = [ E[\det(U^))\9(Y,x) \ ax) = 0]f^ x) (0)dx, 

J B 

where both members are finite. 

Proof. Let B s = {x £ M 3 ;0 < x x < 1,0 < x 2 < f ,0 < x 3 < §} and-B-s = 
B\B 5 . Then 

(7.3) J2 9(Y,x)= £ g(Y,x)+ £ 9(Y,x). 

xgBn^-^o) xeBsr^- 1 ^) aseB-ir^-^o) 

On the set Bs the distribution of £ is nondegenerate, so that Theorem 9.6. 
of [13] can be used. However, that theorem is only valid for a continuous 
function g. Therefore, let {k n } and {h n } be sequences of continuous, mono- 
tone functions such that k n (x) = 1 if x < 0, k n (x) = if x > 1/n, h n (x) = 1 if 
x < — 1/n and h n (x) = if x > 0. It is easy to verify that k n (x) — > l(_ OOi0 i(x) 
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and h n (x) 



L (-oo,0) 



[x) as n — > oo. Because k n (x) and h n (x) are continuous 



for each n the theorem of [13] can be applied with 

g n {Y,x) = k n ( sup W x (s) ) ■ k n (W(xi - x 2 ) - u) ■ k n (W{x 1 + x 3 ) - w) 
\ser x / 



so that 



(7.4) 



x h n (W xx (xi - x 2 )) ■ h n (-W xx (xi + x 3 )), 



E 9n(Y,x) 



E[\det(Ux))\9n(Y,x)\C(x) = 0]f^ x) (0)dx, 



where both members are finite. By dominated convergence arguments as 
n — > oo, (7.4) holds with g n replaced by g. 

On Bs, however, the distribution of £ will degenerate, so the same theorem 
cannot be used in this case. Instead, we show that the assumption of finite 
variance of the number of zeros W x (y) = 0, y £ [0, 1], will force the expected 
value of the first term on the right-hand side in (7.3) to go to zero as 5 
approaches zero. To see that this is the case, first note that g is bounded by 
1 so that 



< E 



E 9(Y,x) 
Lxesn^-Ho) 



E 9(Y,x) 
LzeB-arX-^o) 



< E[JVS(5,)], 



where Nq(Bs) = #{x € Bs;£(x) = 0}. Because W is assumed to be an a.s. 
differentiable Gaussian process, it follows that a.s. there are only finitely 
many solutions to W(y) = in any finite interval. Thus, with probability 
one, {ye[0,l];W(y) = 0} = { yi ,...,yK}, where K = N^([0, 1]) = #{y € 
[0,1]; W(y) = 0}. Using this fact and the stationarity of W, the following 
inequality holds for < 5 < 1 : 



E[Ato)] 



K 



E #{(^3) G [0, 5/2} 2 ; W x { Vl - x 2 ) = W x (vi + x 3 ) = 0} 



li=l 
00 k 



EE E [#{(^3)e(0,<V2] 2 ; 

fc=li=l 

W x ( yi - x 2 ) = W x ( Vl + x 3 ) = 0}}P(K = k) 

00 

E kE[#{(x 2 ,x 3 ) G (0, 5/2} 2 ;W x (-x 2 ) = W x (x 3 ) = 0}}P(K = k) 

k=l 

E[#{(x 2 ,x 3 ) € (0,S/2] 2 ;W x (-x 2 ) = W x (x 3 ) = 0}] • E[N^([0, 1])] 
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<E[#{( s ,t)e[0,i]; 



0<\t-s\<5,W x (s,0) = W x (t,0)=0}]- E[NF([0, 1])]. 

Writing N s = #{(s,t) € [0,1]; < \t - s\ < S,W x (s,0) = W x (t,0) = 0} and 
observing that 

iV ^([0,l])(iV ^([0,l])-l) 

= #{{s,t)e[0,l];s^t,W x (s,0)=W x (t,0) = 0}, 

where N^ x ([0, 1]) = #{y G [0, 1]; W x (y) = 0}, one can by using the formula 
for the second factorial moment, given in [7], deduce that 



E[N 5 ] = / h(s,t)dsdt, 
J{(t,s)£{0,l] 2 fl<\t~s\<6} 

where 

h(s,t) = E[\W xx (s)W xx (t)\ | W x (s) = W x (t) = 0]fw x (.s),W x (t)(Q>Q)- 

Because the assumption of finite variance of the number of zeros N^ x ([0, 1]) 
implies that h(s,t) is integrable over the rectangle [0,1] 2 it follows, by ab- 
solute continuity, that lim^o E[Ng] = 0. Thus 

E 9<y,x) 



= lim / E[|det(6(x))| 5 (r,x) | £(z) = 0]/ c(x) (0) dx, 

so that the final result follows by monotone convergence. □ 

Remark 7.1. In [13] the domain of definition of the function g is C(B,M) x 
B, where C(B,M) is the class of continuous functions on B. However, the 
arguments of that proof holds also in the case when the domain of definition 
of g is C(£>,IR n ) x B, where C(B,W n ) is the class of continuous functions on 
B taking values in M n . The spaces C(B,M) and C(B,M n ) are, in both cases, 
equipped with a uniform convergence norm. 

Remark 7.2. Sufficient conditions for finite variance of the number of 
zeros W x (y) = 0, y G [0, 1], can be found in [7], page 209. 



Using Theorem 7.1, we may state the main theorem for the Palm distri- 
bution of (x2,X3,H2,Hs) in the spatial case. 
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Theorem 7.2. Let the assumptions and notation of Theorem 7.1 hold 
and define 

(7.5) a(x) = W x (x 1 )~W xx (xx - x 2 )~W xx (x 1 + a;3) + l{w a ,( s )<o,Vser iC }- 
Moreover, define a function h s (x,y) = h s (xi,x 2 ,x 3 ,y\, y 2 ) by 

h s (x,y) = E[a(x) \ £(x) = 0,W(xi - x 2 ) = yi, W{x x + x 3 ) = y 2 ] 

x f^{x),W(x 1 -x 2 ),W(x 1 +x 3 ) (0,1/1,2/2) • 

Then the Palm distribution F^ 3 ^ Hi H ^ is given by 

7-ispace / \ 

, . F X2,a> 3 ,H 2 ,H 3 (^ s ^) 

(7.6) 

_ Jo Jq Jo h s (0,X2,x 3 ,y 1 ,y2)dy 2 dy 1 dx 3 dx2 
E[W x (0)~\W(0) = 0]f w{0) (0) 

having the density 

h s (0,r,s,u,w) 



(7-7) fx] P 2 T 3 ,H 2 ,H 3 (r,s,u,w) E[Wx{Q y 1 w{0) = 0]fw(Q){0 y 
where r, s, u > and w < 0. 



Proof. From the definition of the Palm distribution (7.2) and Theorem 
7.1, the numerator in the statement (7.6) of the theorem follows due to sta- 
tionarity in the first coordinate x\. The denominator follows by straightfor- 
ward application of Rice's formula; see the proof of Theorem 6.1. Moreover, 
it follows from (7.6) that _p s P ace is absolutely continuous so that the density 
exists and is given by (7.7). □ 

7.2. Encountered distribution of heights and distances. The correspond- 
ing Palm distribution for overtaking encountered waves will now be consid- 
ered, that is the distribution of distances and heights at times ij > when 
a ship sailing with constant velocity v is overtaken by a wave center. To ob- 
tain a Palm distribution each crossing ij is associated with some distances 
denoted by (x 2 , x 3 )i. In this case (x 2 ,x 3 )i are the distances from vti to the 
closest local maximum and minimum before and after vti, respectively, in 
the wave profile W(x,ti). Following now familiar steps, the following Palm 
distribution can be formed 

F x 2 % 3 ,H 2 ,H 3 (r,s,u,w) 

Efnumber of ti < 1 such that (x 2 ,x 3 ,H 2 ,H 3 ); L < (r,s,u,w)] 
E [number of U < 1] 
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Similarly to the previous case, the Palm distribution can be expressed in 
terms of level crossings of the process W. More precisely, each time U satis- 
fying (x 2 ,x 3 ,H 2 ,H 3 )i < (r,s,u,w) has the following characteristics 

(t,x 2 ,x 3 ) G B, where the set B is defined by 

(CI) 

B = {{t, x 2 , x 3 ) G R 3 ; < t < 1, < x 2 < r, < x 3 < s}, 



(C2) W(vt,t) = W x (vt-x 2 ,t) = W x (vt + x 3 ,t) = 0, 

W xx (vt -X2,t)<0, 

(C3) W xx (vt + x 3 ,t)>0, 

W x (s, t) < 0, Vs G T tx , where Y tx = {s GM.;vt — x 2 < s < vt + x 3 } 
(C4) W(vt-x 2 ,t) <u, W(vt + x 3 ,t)<w, 

(C5) iLwivt, t) = vW x (vt, t) + W t (vt, t) > 0. 



Using these facts and by defining 

(7.8) v (x) = v(t, x 2 ,x 3 ) = (W(vt, t),W x {vt -x 2 ,t), W x (vt + x 3 ,t)) 

the Palm distribution can be written in the following form, suitable for our 
purposes: 

F ™% 3 ,H 2 ,H 3 (r,s,u,w) 

E[#{x 6 B;rj(x) = 0, (C3), (C4) and (C5) satisfied}] 
" E[#{t 6 [0, l];W(vt,t) = 0, (d/dt)W(vt,t) > 0, W x (vt,t) < 0}] ' 
The following theorem is the encountered counterpart of Theorem 7.1. 

Theorem 7.3. Let W :M? — > R 6e a stationary Gaussian process having 
a.s. twice continuously differentiable sample paths. Assume that the spectrum 
of W has a continuous component and that the variance of the number of 
zeros W x (y,0) = 0, y £ [0,1], is finite. From W define a process rj by the 
relation (7.8) and a set B as in (CI). Furthermore, define a vector valued 
process Y = (W, W X ,W XX ), a set Tt x as in (C3) and let g(Y, x) be an indicator 
function defined by 

g(Y,x) = l{w x (s,t)<®ys&r tx }l{w(vt-x2,t)<u}l{w(vt+x 3 ,t)<w} 

X 1 {W xx {vt-X2,t)<V} 1 {W xx (vt+x 3 ,t)>0} 1 {{d/dt)W{vt,t)>Q l }- 
Then, writing Nq{B, g) = G B; n(x) = 0, g(Y, x) = 1}, 

E[NZ(B,g))= [ E[\det( Vx (x))\g(Y,x)\ V (x) = 0]f v{x) (0)dx, 

where both members are finite. 
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The proof of the theorem follows the same lines as the proof of Theorem 
7.1 and is, therefore, omitted. 

We are now ready to state the main theorem for the distribution of {x 2 ,x 3 , 
H 2 ,H 3 ) in the encountered case. The proof is analogous to that of Theorem 
7.2. 

Theorem 7.4. Let the assumptions and notation of Theorem 7.3 hold. 
Define 

(7.10) /3(x)=^(t)+^(^-x 2 ,t)-^ :(; (^ + X3,t) + l { ^ (M) <o,v se r te }, 

where Z(t) = W(vt,t) is the encountered process. 

Moreover, define a function h e (x, y) = h e (xi, x 2 , x 3 , yi, y 2 ) by 

h e (x,y) = E[(3(x) \ V (x) = 0,W(vt - x 2 ,t) = yi ,W(vt + x 3 ,t) = y 2 ] 

x fr)(x),W(vt-x 2 ,t),W(vt+x 3 ,t) (0 5 2/1) i^)- 

Then the Palm distribution F^ C X3 h 2 h z * s 9™ en by 
Fx2% 3 ,H 2 ,H 3 (r,s,u,w) 

(7.11) 

Jo Jo Jo J-00 ^(0, x 2 ,x 3 ,yi,y 2 ) dy 2 dy 1 dx 3 dx 2 
E[Z t (0)+l {Wxm<0} I Z(0) = 0]/ z(0) (0) ' 

having the density 

^ 12 . fx2% 3 ,H 2 ,H 3 ( r i s > u > w ) 

h e (0,r,s,u,w) 



e[^(o)+i { ^ (0) o )<0} I z(o) = o]f m (py 

where r, s, u > and w < 0. 

Earlier it was shown that the relation between spatial and encountered 
distribution of slope could be given a physical interpretation in terms of wave 
velocities. This is also true when it comes to the distribution of distances 
and heights of the local extrema closest to a wave center. To see that this is 
the case, consider the functions a(xi,x 2 ,x 3 ) and (3(t,x 2 ,x 3 ) from Theorems 
7.2 and 7.4, respectively. Due to a relation given in the proof of Corollary 
6.1, namely Z t (0) + l{ l y ;i ,(o,o)<o} = W x (0,0)~ (V — v) + , the functions a and f3 
are related by 

(3(0,x 2 ,x 3 ) = (V - v) + a(0,x 2 ,x 3 ), 

where V = — W t (0, 0)/Wa;(0,0) is the wave velocity at the point (x,t) = (0, 0). 
Moreover, it is easy to see that 7](0,x 2 ,x 3 ) = £(0, x 2 , x 3 ) so that, in fact, 
the functions h e (0,x 2 ,x 3 ,yi,y 2 ) and h s (0,x 2 ,x 3 ,yi,y 2 ) only differ by the 
velocity factor (V — v) + that enters into the expectation in the expression 
for h*. 
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7.3. Evaluation of the densities. So far it has not been mentioned how 
the densities (7.7) and (7.12) can be evaluated. It turns out that for both 
densities the denominators can be computed explicitly. In the first, spatial, 
case, this is straightforward and 



For the denominator of the encountered density (7.12) this is not as easy. 
However, [16] showed that by using certain symmetry properties of the Gaus- 
sian distribution one has 



whenever the spectral moments Aoo> A20, A02 and An are available. 

For the numerators, on the other hand, there are no explicit formulas that 
can be used. Fortunately, there are very efficient numerical routines. In par- 
ticular routines from the Matlab toolbox WAFO (Wave Analysis for Fatigue 
and Oceanography, available gratis at the web page 
http://www.maths.lth.se/matstat/wafo/), custom made for this type of 
calculations can be used; see [6]. However, before any computations can be 
done, the infinite dimensional indicators in the expressions for the densi- 
ties must be approximated by an indicator on a finite set of grid points. In 
both cases, the indicator to be approximated is given by l{w x (u,o)<o,-r<u<s}] 
see (7.5) and (7.10). Therefore, let U = (u±, . . . ,u n ), where — r < u\ < • • • < 
u n < s, be a subdivision of the interval [— r, s] and approximate the in- 
dicator l{W x (u,0)<0,-r<u<s} b y l{W !S (w <) 0)<0 ) i=l ) ...,7i}- In this wa Y> due to tlie 
inequality l{w x (u,o)<o,-r<u<s} < 1 {W x (u z ,o)<o,i=i,...,n}> an u PP er bound of the 
densities is obtained. More information on computational details are given 
in the Appendix. 

8. Examples. The methods presented in this paper can be applied to 
a Gaussian sea having a general spectrum. In the following examples, we 
will evaluate the distributions for the different wave characteristics for a 
longcrested Gaussian sea, that is, a Gaussian sea where the spectrum has 
no angular dependency. To get a longcrested sea let S(u>,6) = S(u>), for one 
single fixed direction 9. In this example 6 = tt, meaning that the waves are 
moving in the direction of the positive rr-axis. The spectrum S(u) will be 
a JONSWAP frequency spectrum; see Figure 2. This is a family of spectra 
fully characterized by the set of parameters (h s ,t p ,^,a a ,o-b), where h s is the 
significant wave height, defined as four times the standard deviation of the 




(7.13) 



E[Z t (0) + l { ^ (0i0)<0} I Z(0) = 0]/ z(0) (0) 
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sea elevation, and t p the peak period. The parameter 7, sometimes called 
the peak enhancement factor, determines the concentration of the spectrum 
around the peak frequency and a a and o\, are spectral width parameters. 
Here, for illustration reasons only, we choose h s = 11.5 m, t p = 12.25 s, 
7 = 1, a a = 0.07 and = 0.09. Because 7 equals one this is also known as 
a Pierson-Moskowitz or Bretschneider spectrum. In this example, we have 
also specified a cut-off frequency u> c , that is, a frequency such that S(u>) = 
for all u > u> c , namely uj c = 1.25 rad/s. 

8.1. Distribution of slope. The distribution of slope observed at centers 
of waves and overtaking encountered waves are given in Theorems 6.1 and 
6.2, respectively. In Figure 3, these distributions are shown, in the encoun- 
tered case for three different ship velocities. Clearly, the encountered distri- 
bution is shifted toward less steep waves and the faster the ship sails, the 
less steep the waves are. A simple explanation to this phenomenon is given 
by the dispersion relation for deterministic waves. Recall that by Corollary 
6.1, the spatial distribution (6.3) is transformed into its encountered version 
(6.6) by the factor (V — v) + . The physical interpretation of this factor is 
that waves that are likely to be slower than the ship will cancel. Because 
the dispersion relation (2.1) implies that steep deterministic waves are slow, 
it thus holds that mainly steep waves are canceled, leading to the shift in 
the distribution — the steep waves are simply not fast enough. However, to 




Frequency [rad/s] 

Fig. 2. The J ON SWAP frequency spectrum used in the examples. This spectrum has sig- 
nificant wave height h s — 11.5 m, peak period t v — 12.25 s and cut-off frequency lo c = 1.25 
rad/s. 
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slope 

Fig. 3. Palm distribution of slope observed at centers of waves in space (solid) and at 
encountered centers of waves overtaking a ship sailing with velocity 7 m/s (dashed), 13 m/s 
(dash-dotted) and 16 m/s (dotted). The computations are done for a J ON SWAP frequency 
spectrum having significant wave height h s — 11.5 m, peak period t v = 12.25 s and cut-off 
frequency ui c = 1.25 rad/s. 

fully understand the relation between slope and velocity, it is too great a 
simplification to look at deterministic waves, although it gives some insight 
into the physics involved in the problem. For a more thorough understand- 
ing, one should study the distribution of the random velocity of the waves 
conditional on the slope. 

8.2. Distribution of waveheight and half-wavelength. Next we consider 
the spatial and encountered joint densities of the waveheight H2 — H% and 
the half- wavelength %2 + £3, that can be obtained from Theorems 7.2 and 
7.4. In Figure 4 these densities are shown, in the encountered case for two 
different ship velocities, namely v = 7, 13 m/s. Clearly, the effect of observing 
the sea surface at waves overtaking the vessel compared to observing it along 
a line at a fixed time, is that the density is shifted toward longer waves in the 
former case. Moreover, the shift increases with increasing ship velocity. As 
for the distribution of slope, this behavior can to some extent be explained 
by the dispersion relation. Recall that [see (4.1)] the velocity of deterministic 
waves satisfying the dispersion relation is proportional to the square root 
of the wavelength. This means that long waves are faster than short ones. 
Because the factor (V — v) + , explaining the difference between the spatial 
and encountered densities, eliminates waves that are slower than the ship, 
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Level curves at: 




50 100 150 



half-wavelength [m] 



Level curves at: , 




} 50 100 150 



half-wavelength [m] 

Fig. 4. Density of half -wavelength X2 + £3 and waveheight H2 — H3 observed at wave 
centers along a line at a fixed time point (top) and at centers of waves overtaking the 
ship (bottom). In the latter case, the velocities of the ship are v = 7 m/s (solid) and 
v — 13 m/s (dash- dotted). The computations are done for a JONS WAP frequency spectrum 
having significant wave height h s = 11.5 m, peak period t p — 12.25 s and cut-off frequency 
lj c = 1.25 rad/s. 
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it thereby eliminates small and moderately high waves in favor of long and 
high waves. 

APPENDIX: COMPUTATIONAL DETAILS 

A.l. Evaluation of densities with WAFO. In this section we discuss the 
evaluation of the densities J s P ace and f cnc from Theorems 7.2 and 7.4, re- 
spectively, using the WAFO-toolbox. 

The density J s P ace in (7.7) can be computed by the routine spec2mmtpdf 
and it computes the desired density given the spectrum of the process 
W(x,t). For the encountered density / cnc given by (7.12), however, there 
is no existing routine in the toolbox, but the function rind that evaluates 
Gaussian multivariate expectations is of great help. 

Before the intensity in the numerator in (7.12) can be evaluated, the 
infinite dimensional indicator must be approximated by an indicator on a 
finite set of grid points of the interval [— r, s]. Let U = (u±, . . . ,u n ), where 
— r < u\ < ■ ■ ■ < u n < s, be a subdivision of the interval [— r, s] and approxi- 
mate the indicator l{ Wx ( u ,o)<o -r<u<s} by l{w x ( Ui ,o)<o,i=i,...,n}- To use rind 
the full distribution, that is, the means and covariances of all variables in- 
volved, must be computed. To simplify notation the variables are collected 
in the following vectors: 

x t = (W X ( U1 , 0), . . . , W x (u n , 0)) = W X (V, 0), 

X d = (W xx (-r,0),W xx (s,0),vW x (0,0) + W t (0,0)), 

X c = (W x {-r, 0),W x (s, 0), W(-r, 0),W(s, 0), W(0, 0)). 

Because W(x,t) is a zero mean stationary process, and due to the fact that 
the means of the (quadratic mean) derivatives of a stationary process are 
zero, the mean of all three vectors above is zero. 

Computing the covariances requires more effort. In general, for a station- 
ary field W(t) with covariance function R(t) the following differentiation 
rule holds 

(^^A^ / y+^s) d c+d W(t) \ +b d a + b + c + d R(r) 

(1.14) ^OV -l-ij 



Let £ be the covariance matrix of (X t ,X^,X c ) and partition it as 



T=t— S 




where it is understood that the lower left part of the matrix is the transpose 
of the upper right part. By repeated use of (1.14) one obtains 

£ tt = (-i^(U T -U,0)), 
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E td = (-R&i-r - U T , 0), -R^(s - U r , 0), 

-^%(-U T ,0)-^ r (-U T ,0)), 

S* c = (-%(-r - U T , 0), -R^(s - U T , 0), -Rs(-r - U T , 0), . . . , 

-R^s-V T ,0),-R^-U T ,0)), 

where — R^(TJ T — U, 0) should be interpreted as the matrix with element 
on row i and column j equal to —R^(ui — Uj,0), and —R^(—r — U T , 0) is 
a column vector with elements —R^(—r — Uj,0), i,j = l,...,n. The other 
submatrices should be interpreted in a similar fashion. The remaining part 
of £ is equal to 



J dc 



/i?««(0,0) R mi (s + r,0) 
i?««(0,0) 



i?^(-r- S ,0) 
-^(-r,0)-^ r (-r,0) 

%(0,0) 
%(-r-s,0) 



«ii^(r,0) + .% T (r,0) 
t; J R eee (-s,0)+ J R eer (-s,0) 
-u 2 .%(0,0) - 2^ r (0,0) - i? TT (0,0) 

i2e«(* + r,0) 


-vR a (s,0) - Rt T (s,0) 



ui? e (-r,0) -i? r (-r,0) -vi? c (s,0)-i? T (s,0) 





-i?«(0,0) -i^(s + r,0) 
-7^(0,0) 



-i? 5 (-r- S ,0) 
i?(0,0) 



-^(s + r.O) 


7?(s + r,0) 
#(0,0) 



V 



J2«(r,0) 
%(-s,0) 


-^(r-,0) \ 
-7? e (- S ,0) 
E(r,0) 
i?(-s,0) 

i?(0,0) y 



The covariance function and its derivatives can be computed by the WAFO- 
function spec2cov, that from a given spectrum computes the covariance 
function and its derivatives up to the fourth order. 



A. 2. Computational issues. When the conditional expectation in the nu- 
merator of (7.12) is computed numerically new difficulties arise. If the co- 
variance matrix is badly scaled, densities can become singular in the nu- 
merical computations. For example, it might happen that variances become 
negative even though they should be small and positive. A way to make 
the algorithms more stable is to rescale the covariance matrix such that its 
elements are of comparable sizes. 

For the spectra we will consider the average period is about 10 seconds 
and the value of the average wave length is about 10 times as large, but 
measured in meters. Because the sea elevation takes values in the same range 
in both time and space, the derivatives will be of very different magnitudes. 
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To get a better scaled covariance matrix the time- and space-coordinates 
are, therefore, transformed according to 





t 

Making this change of variables the spectral density is transformed to 

s{{v,e) = s{oj,( 




Note that the angle 9 remains unaffected. In this way, the spectral moments 
Aoo> A20 and A02 of S(u),6) all become equal to one; see [17] for a proof. This 
means that the average number of waves per transformed meter is equal 
to 1/2tt as well as the number of waves per transformed second. Also, the 
variance of the sea elevation is transformed to one. Note that it is common 
practice in wave data analysis to normalize the data to have unit variance. 
However, in this case this is not enough in order to avoid a badly scaled 
covariance matrix because the problems arise due to different spatial and 
temporal scales. When the coordinates are transformed, the velocity of the 
ship changes to 




To relate the spectral moments to the elements in the covariance matrix note 
that according to (2.2) and (2.3) it holds that A 00 = R{0, 0), A 20 = -%(0, 0) 
and A02 = — i2 rr (0, 0). Because many of the elements on the diagonal of 
£ are of this form, we can conclude that the covariance matrix obtained 
after the coordinate transformation is better scaled than the original one. 
After the density is computed with the transformed coordinates, it is then 
easily transformed back to the true units seconds and meters. These kind of 
transformations are implemented in the function wnormspec in WAFO. 
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